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ABSTRACT 

I discuss the effect of non-radial motions on the small-scale pairwise peculiar 
velocity dispersions of galaxies (PVD) in a CDM model. I calculate the PVD for the 
SCDM model by means of the refined cosmic virial theorem (CVT) (Suto & Jing 1997 
(hereafter SJ97)) and taking account of non-radial motions by means of Del Popolo & 
Gambera (1998) (hereafter DG98) model. I compare the results of the present model 
with the data from Davis & Peebles (1983), the IRAS value at l/i _1 Mpc of Fisher 
et al. (1993) and Marzke et al. (1995). I show that while the SCDM model disagrees 
with the observed values, as pointed out by several authors (Peebles 1976, 1980; Davis 
& Peebles 1983 (hereafter DP83); Mo et. al 1993; Fisher et al. 1994b; SJ97; Jing et 
al. 1998 (hereafter J98)), taking account of non-radial motions produce smaller values 
for the PVD. At r < lft, _1 Mpc the result is in agreement with Bartlett & Blanchard 
(1996) (hereafter BB96). At the light of this last paper, the result may be also read as 
a strong dependence of the CVT prediction on the model chosen to describe the mass 
distribution around galaxies, suggesting that the CVT cannot be taken as a direct 
evidence of a low density universe. Similarly to what shown in Del Popolo & Gambera 
(1999, 2000) (hereafter DG99, DGOO), Del Popolo et al. (1999) (hereafter D99), the 
agreement of our model to the observational data is due to a scale dependent bias 
induced by the presence of non-radial motions. Since the assumptions on which CVT 
is based have been questioned by several authors (BB96; SJ97), I also calculated the 
PVD using the redshift distortion in the redshift-space correlation function, £,z(r p ,n), 
and I compared it with the PVD measured from the Las Campanas Redshift Survey 
by J98. The result confirms that non-radial motions influence the PVD making them 
better agree with observed data. 
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1 INTRODUCTION 

The pairwise velocity dispersion of galaxies (PVD) is an important quantity that gives information on the structure and 
clustering of the universe. Peculiar velocities, originating from the action of gravitational fields, are a probe of the gravitational 
potentials produced by luminous and dark matter. 

Several approaches have been developed to determine the PVD. In this paper, I'll be concerned with two of them, the 
first one based on the redshift-space distortions and the second one on the Cosmic Virial Theorem (hereafter CVT). 

The first method relies on the fact that peculiar motions do not affect tangential distances, but only the radial ones. If 
we observe the three-dimensional distribution of galaxies in redshift-space, we can see that it is distorted with respect to that 
seen in real space: dense clusters appear elongated along the line of sight ('fingers of God' effect) and the correlation functions 
of galaxies are also distorted when they are viewed in redshift-space. The important point is that the redshift-space distortion 
effects can be used to get useful information on important quantities like the PVD. This can be done using the anisotropics 
created by the peculiar motions in the redshift space correlation function. From a redshift survey of galaxies, it is possible to 
measure the correlation function in redshift-space, £(s). One can calculate £(s) as function of the component of the separation 
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vector in the direction parallel (n) and perpendicular (r p ) to the line of sight. Plotting the correlation function £z(r p ,n) in 
the 7r-r p plane, redshift distortions give rise to: 

1) a stretching of the contours of £z(r p ,ir) along n on small scales (< a few Mpc) because of non-linear pairwise velocities; 

2) a compression along the line of sight, n, on larger scales, because of bulk motions. 

DP83 used the small scale distortion to get the PVD of galaxies, modelling the elongation of fz(r Pl 7r) along n as a convolution 
of the real-space correlation function with the distribution function /(V12), being V12 the relative velocity of the galaxy pair 
in the direction of the line of sight. The quoted method, introduced by Davis et al. (1978), Peebles (1980), DP83 and Bean et 
al. (1983) has sometimes given discrepant values: DP83, using data from the CfAl, obtained a value of o\2 = 340 ± 40km/s, 
at r = l/i _1 Mpc; Fisher et al. (1994a) using the IRAS 1.2 Jy survey, obtained a value of CT12 = SlTl^km/s in agreement with 
that of DP83. Larger values was obtained by Marzke et al. (1995) using the CfA2 and the SSRS2 (Southern Sky Redshift 
Survey) redshift surveys, (<7i2 = 540 ± 180km/s), by Lin et al. (1995) (0-12 = 452±60km/s) using the Las Campanas Redshift 
Survey and by J98 (<ti 2 = 570 ± 80km/s). Moreover, some authors (Mo et al. 1993; Marzke et al. 1995; Guzzo et al. 1997) 
have pointed out that the value of the PVD is very sensitive to the presence or absence of rich clusters in a sample. Other 
authors emphasized that being o\2 a pair-weighted statistic it is heavily weighted by the densest regions of a sample (Strauss 
1997). In two recent papers, Sheth et al. (2000) and Jing et al. (2001) have shown that (T12 varies widely not only because 
of the previous effect, but also because of the morphological type of the galaxies entering the catalog. In particular, Sheth et 
al. (2000), showed that redshift space distorsions should affect red galaxies more strongly than blue. This result can partly 
explain the difference between the IRAS catalog o\2 and the CT12 computed from the more fair sample contained in the Las 
Campanas Redshift Survey (LCRS). However the most recent papers agree to values of the PVD in the range 400-600 km/s 
(J98). This final range of convergence of cri2, although still large, comes because of the increasingly large volumes surveyed 
in the most recent work. 

The (T12 statistics has been used as a discriminator of models. DP83's result compared with the SCDM model showed 
that the model predicts values of PVD larger than the observed ones and successive studies confirmed that the SCDM model 
disagrees with the observed values by a factor 3 — 6 (Mo et al. 1993; Fisher et al. 1994b; Marzke et al. 1995). A new idea, 
that of biased galaxy formation (Kaiser 1984; Bardeen et al. 1986 (hereafter BBKS)), was needed to reduce the discrepancy 
between the SCDM model prediction and observations. Some recent studies have however shown that for the SCDM model 
the shape of the predicted PVD is different from that observed, since the SCDM model predicts a correlation function steeper 
than the observed one (J98). In any case, introducing a constant bias does not resolve the problem, and a scale-dependent 
bias is necessary to reduce the quoted discrepancy (J98). 

Another tool that has long been used to determine o\i is the Cosmic Virial Theorem (CVT) (Peebles 1976) relying on 
the assumption of hydrodynamic equilibrium of galaxies on small scales. Several authors (Fisher et al. 1994b; BB96; SJ97) 
have pointed out some limitations of the theorem. BB96 discussed how the mass distribution around a typical pair of galaxies 
affects the CVT predictions, while SJ97 improved the theorem taking into account the effect of the finite size of galaxies. 
These last authors showed that the CVT in its initial formulation overestimates the small-scale PVD. 

In this paper, I'll study the effects of a non-local bias, induced by non-radial motions (see DG98, DG99, DG00; D99), on 
the PVD by means of both the CVT and the distortions in £ z (r p ,7r). 

The plan of the paper is the following: in section 2, I calculate the PVD using the CVT in the original form introduced 
by Peebles (1976) and the improved version of SJ97 and taking account of non-radial motions. The theoretical results are 
compared with observational data. In section 3, I calculate the PVD by means of the distortions in the redshift-space 
correlation function and compare the results with observations. Section 4 is devoted to conclusions. 



2 PVD AND THE CVT 

2.1 The CVT and the finite size of galaxies 

The CVT introduced by Peebles (1976) is fundamentally an equation of hydrostatic equilibrium obtained assuming statis- 
tical equilibrium between galaxy pairs, on small scales. The formula describing this theorem, and giving the relative (one- 
dimensional) peculiar velocity dispersion as function of their separation r, can be obtained by means of the second BBGKY 
hierarchy equation. It can be expressed as: 



where is the mean density of the universe, £ p and C, p are the two and three-point correlation function of mass, respectively. 
Assuming the forms for the two and three-point correlation function of galaxies given by Groth & Peebles (1977) and DP83, 
namely: 
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Figure 1. The threshold 8 C as a function of the mass M, through the peak height u, taking account of non-radial motions with the model 
of this paper (dashed line) is compared with the result of Sheth et al. (1999) (solid line), obtained using an ellipsoidal collapse model. 



and 

C t (n,r a ,^) = g K e (ri)&(r 3 )+6 s (r3)&(r3)+&(r 8 )&(ri)] (3) 

where Q g = 0.7 ± 0.21 (see BB96), r = (5.4 ± 0.3)/i _1 Mpc, 7 = 1.77 ± 0.04, and assuming that the two and three-point 
correlation functions of mass, £ p and £ p follow the same scaling and satisfy the relations: 

Ur) = b 2 Ur) (4) 
and 

Cg(n, r 2 ,r 3 ) = Q s b A C,p(r 1 ,r 2 ,r 3 )/Qp (5) 
being b the bias parameter, the CVT gives (Peebles 1976; Suto 1993): 



K 2 M)- = 1460W?#,/§M5.4^ U-*) \^-\ Ws (6) 



2— y 



1.3b 2 V 33.2 ' ybAh-iMpc J ^l/i^Mpc 

where -^0(7) is a function given in Peebles (1980) (equation (75.12)) and fio is the density parameter. Numerical and analytical 
studies in non-linear gravitational clustering seem to indicate that Q p can be approximated as a constant in the range 0.5 — 2 
(Suto 1996), which depends very weakly on the underlying cosmological model. We shall use the value 0.6, that gives the same 
predictions of BB96 uncorrected CVT (see their Fig. 1), which shall be used for comparisons with our paper (note that also 
a low value of this parameter gives prediction for the pairwise-peculiar velocity dispersions much larger than observations). 

Equation (^|) was improved by SJ97 by taking account of the fact that galaxies are not point-like, using perturbation 
analysis and numerical integration. They changed the CVT incorporating a non-zero core radius, r c , and softening the 
gravitational force by means of a Plummer law with softening radius r s obtaining the final result: 



Mr)) = 1460 ^/ w 5A~^ 



5.4h -1 Mpc 

I(r c /r, r s /r;7) 



33.2 V lft _1 Mpc 
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Figure 2. The bias factor b(v) as function of i/ 2 . The solid line represents the spherical collapse prediction of Mo & White (1996), the 
dotted line the prediction for b obtained from the model of this paper and the dashed line the ellipsoidal collapse prediction of Sheth et 
al. (1999). 

where I(r c /r, r s /r; 7) is given in SJ97 (equation (10)). This result simply tells that by replacing /o(t) with I(r c /r, r s /r; 7) the 
finite size of the galaxies is automatically taken into account. Moreover they showed, using COBE normalized CDM models, 
that the quoted improvement of the CVT has as final result a reduction of the small-scale velocity dispersion of galaxies, but 
not enough to be in agreement with observations. 

The final aim of this section is to determine the PVD, using the modified CVT and, moreover, taking into account 
non-radial motions. 

2.2 The threshold S c , the selection function and bias 

The goal quoted at the end of the previous subsection can be accomplished remembering that non-radial motions induce a 
non-local scale dependence of the bias parameter, b (see D99)) and taking into account this effect in equation (0). 
In this subsection, we are going to show how to calculate the dependence of b on the peak height v = ^1 j which is proportional, 
through the overdensity S = p ~^ b and a, the rms value of 8, to the halo mass. The b(v) dependence obtained shall be used 
in the next subsection to analyse the effects of non-radial motions on the PVD. 



As shown by DG98, DG99, DG00, if non-radial motions are taken into account, the threshold S c is not a constant but is 
function of mass, M, (DG98, DG99, DG00): 



where 8 co — 1.68 is the critical threshold for a spherical model, n is the initial radius, L the angular momentum, Ho and fio 
the Hubble constant and the density parameter at the current epoch, respectively, a the expansion parameter and 8 the mean 
fractional density excess inside a shell of given radius and mass M. The mass dependence of the threshold parameter, S c (v), 
and the total specific angular momentum, h(r, v) — L(r, v)/M, acquired during expansion, were obtained in the same way as 
described in D99 (see also Fig. 1 of the quoted paper). 

The result of the calculation is shown in Fig. 1, where I plot S c (v) obtained by means of the model of this paper together with 
that obtained by Sheth et al. (1999) using an ellipsoidal collapse model. The dashed line represents 8 c (v) obtained with the 
present model, while the solid line that of Sheth et al. (1999). Both models show that the threshold for collapse decreases with 
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mass. In other words, this means that, in order to form structure, more massive peaks must cross a lower threshold, 8 c (v), with 
respect to under-dense ones. At the same time, since the probability to find high peaks is larger in more dense regions, this 
means that, statistically, in order to form structure, peaks in more dense regions may have a lower value of the threshold, 8 c (y), 
with respect to those of under-dense regions. This is due to the fact that less massive objects are more influenced by external 
tides, and consequently they must be more overdense to collapse by a given time. In fact, the angular momentum acquired 
by a shell centred on a peak in the CDM density distribution is anti-correlated with density: high-density peaks acquire less 
angular momentum than low-density peaks (Hoffman 1986; Ryden 1988). A larger amount of angular momentum acquired by 
low-density peaks (with respect to the high-density ones) implies that these peaks can more easily resist gravitational collapse 
and consequently it is more difficult for them to form structure. This is in agreement with Audit et al. (1997), Peebles (1990) 
and DP98, which pointed out that the gravitational collapse is slowed down by the effect of the shear rather than fastened 
by it (as substained by other authors). Therefore, on small scales, where the shear is statistically greater, structures need, on 
average, a higher density contrast to collapse. 

This results in a tendency for less dense regions to accrete less mass, with respect to a classical spherical model, inducing 
a biasing of over-dense regions towards higher mass. Moreover the scale dependence of the threshold, 8 c (v), implies also a scale 
dependence of the bias parameter, b, since the two parameters are connected (see Borgani 1990; Mo & White 1996; DG98, 
D99). In fact, according to the biased theory of galaxy formation, observable objects of mass ~ M arise from fluctuations of 
the density field, filtered on a scale Rf, rising over a global threshold, 8 > S c = vta, where a is the rms value of 8 and Vt is the 
threshold height. The number density of objects, n p k, that forms from peaks of density of height v can be written, following 
BBKS, in the form: 

P OO 

n pk = / t(—)N pk {v)dv (9) 
Jo Ut 

where t(-j^) is the threshold function, and N^dv the differential number density of peaks (see BBKS - equation (4.3)). The 
threshold level i/t is defined so that the probability of a peak becoming an observable object is 1/2 when v = vt. In the sharp 
threshold case the selection function, is a Heaviside function t(~) = &( v — The threshold function is connected to the 
bias coefficient of a class of objects by (BBKS): 

b(Ri) = M + i (10) 
where (y) is: 

70 



1 -7 2 



t(—)N pk (v)dv (11) 



while, 7 and i? are given in BBKS (respectively equation (4.6a); equation (6.14)). 

The threshold (or selection) function can be obtained following the arguments given in Colafrancesco et al. (1995) and DG98. 
In this last paper the selection function is defined as: 

/•OO 

t(y)= / p\S,(5(rMt,u)),aj(rMt,u)]dS (12) 



where S is the mean fractional density excess inside a given radius, as measured at the current epoch, assuming linear growth, 
< 8 > its average value (see Ryden & Gunn 1987) and oj its dispersion given in Lilje & Lahav (1991), and where the function 

gives the probability that the peak overdensity, 8 is different from the average, in a Gaussian density field. The selection 
function depends on v through the dependence of S(r) on v. As displayed, the integrand is evaluated at a radius ryit which is 
the typical radius of the object that we are selecting. Moreover, the selection function t{v) depends on the critical overdensity 
threshold for the collapse, 8 c (y). Given 8 c {v) and chosen a spectrum, the selection function is immediately obtained through 
equation ( |l2] ) and equation (|l3|). As shown in DG98 (see their Fig. 6), the selection function, as expected, differs from an 
Heaviside function (sharp threshold). The shape of the selection function depends on the values of the filtering length Rf and 
on non-radial motions. The value of v at which the selection function t(y) reaches the value 1 (t(u) ~ 1) increases for growing 
values of the filtering radius, Rf. This is due to the smoothing effect of the filtering process. The effect of non-radial motions 
is, firstly, that of shifting t(u) towards higher values of v, and, secondly, that of making it steeper. 

This means that while in a 8 threshold scheme, fluctuations below 8 c (u) have zero probability to develop an observable 
object, and fluctuations above 8 c (y) have zero probability not to develop an object, the situation is totally different when the 
threshold is not a sharp step function. In this case objects can also be formed from fluctuations below 8 c (y) and there is a 
non-zero probability for fluctuations above 8 C to be sterile. This result, in agreement with Borgani (1990), with the 'fuzzy' 



© RAS, MNRAS 000, §-?? 



6 A. Del Popolo 



threshold approach of Audit et al. (1997) and Sheth et al. (1999), is fundamentally connected to non-sphericity effects present 
during the gravitational growing process. 

I want to remark that the choice made to calculate 6, instead of the well known Mo & White (1996) model: 

6=1+ (14) 

is motivated by the fact that I am interested in studying the PVD on scales < 10h~ Mpc, while in their model, Mo & White 
(1996) obtained analytic results only in the limit of large separations (Sheth & Lemson 1999, Taruya & Suto 2000). Although 
not directly necessary to the development of the remaining part of the paper, I also calculate the large scale bias factor, b, in 
order to test the model described and to compare it with the value found by Sheth et al. (1999), as previously done with the 
threshold, <5 C . 

In Fig. 2, I plot the bias parameter, 6, as a function of the peak height u, which is proportional to the halo mass. The 
solid line shows the spherical collapse prediction of Mo & White (1996), the dotted line the prediction for b obtained from 
our model, and the dashed line the ellipsoidal collapse prediction of Sheth et al. (1999). As shown in the figure, the effect of 
non-radial motions is to change the dependence of 6 on v in good agreement with Sheth et al. (1999). From Fig. 2 it is evident 
that at the low mass end the bias relation has an upturn, meaning that less massive haloes are more strongly clustered than 
the prediction of the spherical collapse model of Mo & White (1996) and in agreement with N-body simulations (Jing 1998; 
Sheth & Lemson 1999). 

2.3 PVD, CVT and non-radial motions 

At this point I am ready to calculate the PVD, using the CVT in its original formulation (Peebles 1976) and the improved 
version of SJ97, and to study the effect of non-radial motions on these results. 

Fig. 3 plots the results of the model of this paper and the comparison with the DP83 data for values of the parameter 
F = 0.1, 1, 1.5 (open squares), the IRAS value at l/i _1 Mpc (filled square) (Fisher et al. 1993), the DP83 interpretation of the 
Turner (1976) galaxy pair catalog (solid hexagons) and the Marzke et al. (1995) data at l/i _1 Mpc (dashed errorbar). Q Here, 
the PVD, is derived using a SCDM model (Qo — 1, h — 1/2) filtered on galactic scales and normalized to Qcobe = 17/iK. 
The dotted line represents the PVD obtained from the uncorrected version of the CVT. 

As previously found by several authors (DP83; Mo et al. 1993; Fisher et al. 1994b; Marzke et al. 1995), the SCDM model 
overestimates the PVD. The dashed line represents the PVD obtained from the corrected version of the CVT (SJ97), having 
assumed r c = r B = 10/i _1 kpc. 

In agreement with SJ97, taking into account the finite size of galaxies reduces the small-scale velocity dispersion of 
galaxies, but not enough to be in agreement with observational data. The solid line represents the PVD, taking account of 
the non-radial motions, while the long-dashed line takes also into account the finite size of galaxies (again r c = r s = 10/i -1 
Kpc). The finite size effect and non-radial motions both produce a decrease in the PVD. This is due to the fact that the 
finite size effect suppresses the effective gravitational force between pairs and changes the two-point correlation function on 
small scales. Non-radial motions have also the effect of changing the two-point correlation function (DG99; D99), the mass 
distribution and the density profile of the galactic halos (White & Zaritsky 1992; DG98). The change of mass distribution 
produces, in agreement with BB96, a change in the PVD predicted by the CVT, namely smaller values of the small-scale 
pairwise peculiar velocity dispersions of galaxies. Differently from BB96 in the model of this paper, I do not observe the rapid 
increase of the PVD beyond r ~ l/i _1 Mpc with a consequent better agreement of the model with data. The difference is due 
to the increase of the bias parameter going from small to larger scales. Aside from the model introduced to calculate b, the 
behavior of the bias parameter, and then the PVD, can be qualitatively explained as follows: the papers of Sheth (1996) and 
Diaferio & Geller (1996) and Mo et al. (1997) suggest that, while on very small separations, pairs come from both small and 
large haloes, at larger separations, pairs come mainly from larger and more massive haloes. Since massive haloes have larger 
b values than less massive ones, this means that the effective bias is lower on small scales than on larger ones. 

The CVT has been traditionally applied as an indicator of the cosmological density parameter, Qo , and since its first use 
(DP83), it was considered as a strong indicator of an open universe. Excluding the behaviour of the PVD at r > l/i _1 Mpc, 
where, as remarked by BB96, the stable clustering hypothesis could be questioned, the result of this paper is in agreement 
with that of BB96, namely: the PVD should not be taken as a direct evidence of a low density universe. At the light of Fig. 
3, I want to add that this conclusion should be more strict if we could rely on the CVT for r > l/t -1 Mpc. At the same 
time, as remarked by SJ97, the usefulness of the CVT to put constraint on fio is strictly connected with the quality of the 

* We recall that the parameter F was introduced in the infall model used by DP83 (see DP83 equation (22)- (23)) in order to test the 
sensitivity of the velocity dispersion estimates to the streaming correction. According to the equation (23) in DP83, a small value of F 
implies that galaxies on all scales expand with the Hubble flow, while a value of F = 1, implies, that galaxies having separation of 5Mpc 
are expanding at 1/2 the Hubble rate. 
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Figure 3. The PVD as a function of separation r. The PVD, is derived using a SCDM model (Qq = 1, h = 1/2), whose transfer function 
is given in next section, equation (pr^), filtered on galactic scales and normalized to Qcobe = 17/i K. The dotted line represents the PVD 
obtained from the uncorrected version of the CVT, the dashed line represents the PVD obtained from the corrected version of the CVT 
(SJ97) with r c = r a = 10h~ ^pc, the solid line the PVD taking into account the non-radial motions while the long-dashed line takes 
also into account the finite size of galaxies (r c = r s = 10/i -1 kpc). The model is compared with DP83's data for values of the parameter 
F = 0.1,1,1.5 (open squares), the IRAS value at lh~ Mpc (filled square) (Fisher et al. 1993), the DP83 interpretation of the Turner 
(1976) galaxy pair catalog (solid hexagons) and the Marzke et al. (1995) data (dashed errorbar). 



observational data. In fact, as reported in the introduction, some authors (Mo et al. 1993; Marzke et al. 1995; Guzzo et al. 
1997) have pointed out that the value of the PVD is very sensitive to the presence or absence of rich clusters in a sample. 
There is a large variation in < vf 2 > in different surveys and even for the same survey analyzed in different ways. Moreover, 
<7i2 varies widely also because of the morphological type of the galaxies entering the catalog (Sheth et al. 2000). As shown by 
Mo et al. (1997), the numbers of pairs with small separations is dominated by haloes having M ~ M* while on small scales 
massive pairs, M > M,, strongly affect the PVD. Since the number density, n, of this massive haloes is small, a large sample 
is needed to have a fair sample of n or in other terms the small scale PVD can be fairly sampled only by means of samples 
containing many rich clusters. 



3 PVD FROM DISTORTIONS IN £z(-Rp,7r) 

Since the combination of statistical and systematic uncertainties entering the CVT have lead several authors to conclude that 
it has notheworthy problems as an estimator of cosmological parameters (Fisher et al. 1994b), in order to be on the safe side, 
I also use a more reliable and more widely used method to estimate cosmological parameters, and in particular the PVD. 

As previously reported in the introduction, the three-dimensional distribution of galaxies in redshift-space appears dis- 
torted with respect to the same distribution in real space. The correlation function measured in redshift space, £(s), is different 
from the real space counterpart, £(r), because of two effects. These effects are seen in redshift-space, plotting the correlation 
function in terms of two variables, the separations parallel (n) and perpendicular (r p ) to the line of sight. Given a pair of 
galaxies with redshifts corresponding to velocities vi and V2, the separation in redshift space is given by: 

S = Vi - v 2 (15) 

while the observer's line of sight is: 
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l=i(v 1+ v 2 ) (16) 
The separations parallel and perpendicular to the line of sight are respectively: 

* = I (17) 
and 

r p = ss - TV 2 (18) 

The PVD can be obtained by modelling the redshift distortion of £ z (r* p , n), as follows: 

The statistics £ z (r p ,7r) is a convolution of the real-space correlation function, £(r), with the distribution function of the 

relative velocity along the line of sight, f(vvij- Following, for example, Fisher et al. 1994b, Jing & Borner (1998) (hereafter 
JB98), we have: 



l + £.(r P ,ir) = / /(«i 2 ) 



J f(vi2) [l + £ (y/r* + (ir-v 12 /H y 



dv 12 (19) 



To obtain the PVD from £ z (r p ,7r) there are at least three methods (Peebles 1980, section 76). The first one, which I am 
going to use in this paper and which is the more diffused in papers that try to recover the PVD from survey of galaxy redshift 
(Davis & Peebles 1983; Mo et al. 1993; J98; JB98), can be summarized in the following steps: 

1) Estimate of the redshift-space two-point correlation function £ z (r p , n). If one has observational data, coming from a redshift 
survey, it is necessary to choose an estimator for £ z (r p ,7r) (see for example J98). In our case, £ z (r p ,7r) is to be calculated 
theoretically from the power spectrum, as I am going to do. 

2) Estimate of the projected two-point correlation function w(r p ): 



*»(r p ) = 2 Ur P ,K)dir = 2 £(y/r* +y*)dy (20) 



3) Calculation of the two-point correlation function £(r). 

4) Assumption of a functional form for /(V12), an infall model for V12 to, finally, solve equation ( |l9|) for <7i2. 

For seek of completeness, I shortly describe the other two methods. The second approach follows from Eq. |l9[ After 
determining £(r) from: 



C(rp,7r)= / Z(y/r$ + v 2 )dy (21) 
Jo 

the PVD can be obtained by the equation: 

/ / Ur P ,*)dK (22) 



< «12 >= 3Hq 



(see also Peebles 1979). This method does not require the knowledge of f(vi2) but, as noted by Davis & Peebles (1983), the 
integral depends sensitively on £ at \tt\ > lOOOkm/s, where, if one uses data from redshift surveys, £ z (r p ,7r) is poorly known. 
The third method also deriving from Eq. [it] is: 

/•oo poo 

<v 2 2 >=3H 2 / dirdr p (-ir 2 - r p )£ z (r p , it)/ / £z(r p , Tv)dndr p (23) 
Jo Jo 

which requires only £ z (r p ,7r) but also this approach, as the previous one, is not free from drawbacks (see Peebles 1980, Sect. 

76). 

Coming back to the first method, we can start from the step (1). 

In order to get £z(r p ,7r), I Fourier transform the redshift-space power spectrum: 

£(s,A*si) S = 7^3 y"d 3 kP S (fc,^ kl )exp(-ifcs Mks ) (24) 

(see Cole et al. 1994, equation (Bl)), where, as previously quoted, 1 denotes the direction of the line of sight, /iks and 
Mki — ( fc 2 +fc 11 )i/2 the cosine between the vectors s, k and k, 1, respectively (k± and k\\ axe the components of the wavevector 

perpendicular and parallel to the line of sight). £ z (r p ,7r) can be obtained from equation (JmJ) by changing variables according 
to equation 

Before going on, I want to recall that there are three effects that cause a departure of the observed clustering properties 
of galaxies from the linear spectrum: 

a) Redshift space effects. We extract three dimensional clustering information from redshift surveys; in these surveys, the 
galaxy radial coordinates are distorted by peculiar velocities. 
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b) Nonlinear evolution, producing, on small scales, a departure of the mass spectrum from its initial form. 

c) Bias. Different classes of objects trace mass in different ways leading to difficulties in connecting theory and observations. 

Firstly, (relatively to point a), in order to calculate the redshift space power spectrum, distorted by the peculiar velocity 
field, I use the expression given by Peacock & Dodds (1994) and Cole et al. (1995): 

P S (fcx,fc||) = P R (k) [l + /?Mki] 2 (i + fcV 2 2Mki/2) -2 (25) 



where k = + k^_, P R (k) is the real-space power spectrum, and /3 = where the perturbation parameter /(f2) is 

defined in Peebles 1980 (section 14). The real-space power spectrum that I adopt is P(k) = AkT 2 (k) with the transfer 
function T(k) given in BBKS (equation (G3)): 

T(k) = [ln(1 2 + 3 y 4g)] • [1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71) 4 p 1/4 (26) 

where A is the normalizing constant and q — n ^fil^' pc -i ■ Here 8 = p cr /(1.68p 7 ) represents the ratio of the energy density in 
relativistic particles to that in photons (9 = 1 corresponds to photons and three flavors of relativistic neutrinos). The power 
spectrum was COBE normalized using the cosmic microwave anisotropy quadrupole Qrms-ps = 17/iK that corresponds to 
as = 0.95 ± 0.2 (Smoot et al. 1992; Liddle & Lyth 1993). 

Secondly, (relatively to point b), since I have to calculate the PVD in the range 0.01 < r p < 10h -1 Mpc I should take 
into account the non-linear evolution of the power spectrum. 

To this aim, 1 shall use the fitting formula of Peacock & Dodds (1996) (see also Hamilton et al. (1991), Peacock & Dodds 
(1994), Jain et al. (1995), Mo et al. (1997)), relating the evolved power variance, A E = ^^k s P(k)E, to the initial density 
spectrum: 

A|(fc E ) = F([A 2 (k L )]) (27) 

where L stands for 'linear' E for 'evolved', &l = [1 + A E (fcE)] _1 ^ 3 fcE and the functional form of F is given in Peacock & Dodds 
(1996). This correction was done to the real-space spectrum in equation ( p5j ) which shall be introduced in equation ([mJ) to 
obtain £ z (r p , n). 

After £ z (r p ,7r) has been calculated, we can go to step (2) and (3). 

Following the usual technique used to recover the real-space correlation function from (, z (r p ,iv), I define a projected 
function w(r p ), unaffected by redshift distortions: 



w(r p )=2 Ur P ,K)dv = 2 t(^r 2 +y 2 )dy (28) 
Jo Jo 

(DP83). Then £(r) can be obtained, following DP83, either solving w(r p ) for £(r), or also fitting w(r p ) to a power law model 
for £(r). Following the first alternative £(r) is given by evaluating the Abel integral: 

1 f°° 

€W = -^y dr p w(r p )(r 2 p ~r 2 r 1/2 (29) 

(DP83) (an alternative method to calculate £(r), even if more complex, and compatible with that previously presented, as 
remarked by Peebles (1983) and J98, should be that presented in D99, and based on the use of the Limber's equation. As 
shown in that case, the effect of non-radial motions is that of producing enough additional clustering to fit the £(r) of the 
APM galaxy survey). 

Finally, (step 4), I need a functional form for /(V12). A good approximation to this function is: 
}{v 12 ) = -= exp ( ! '-) (30) 

V2(T12 °"12 

(Fisher et al. 1994b; J98), where 1J12 and ai2 are the mean and the dispersion of the 1-D pairwise peculiar velocities. The infall, 
V12, is difficult to model because it is correlated with ai2, and it is scale dependent. Following J98, 1 assume a self-similar 
infall model for 1J12: 

^ 12 M = 1 ^2 ( 31 ) 

where r* = 5h~ Mpc and y is the radial separation in real space. This infall model is usually assumed, because it gives a 
good approximation to the infall pattern seen in CDM models and, as shown, in JB98, J98 the PVD reconstructed using this 
model is in good agreement with the 3-D velocities in the simulations. 

As stressed by J98, the method described to determine the PVD is approximated for several reasons: 
1) the functional forms of /(V12) and Vi2{r) are approximations; 
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Figure 4. The PVD calculated from the redshift-space distortion of £ z (rp,7r). The solid line is the PVD taking account of non-radial 
motions, the dashed line the PVD obtained from the SCDM model that does not take account of the non-radial motions. The data are 
the PVD measured from the Las Campanas Redshift Survey by J98. 



2) the reconstruction of the PVD from redshift distortion gives, as a result, only an average of the real PVD along the line of 
sight. 

In any case, the value of (T12, obtained with the quoted technique, is within 20% of the true PVD (JB98). 

At this stage, I have all the quantities required by equation ([[9]) to obtain £z(r PI n). I have obtained the PVD by using 
equation ( |l9| ) and fitting eri2 to £ z (r p ,7r) (see Peebles 1980; DP83; Mo et al. 1993; J98). The calculation was repeated two 
times, the first with S c = 1.686 and the second assuming equation (]8]). 

The result of the model described in this section is shown in Fig. 4. Here the solid line represents the PVD, taking into 
account non-radial motions, the dashed line the PVD obtained from the SCDM model that does not take into account the 
non-radial motions. The theoretical results are compared with the PVD measured from the Las Campanas Redshift Survey by 
J98. As shown in Fig. 4, the PVD predicted by the SCDM model that does not take into account non-radial motions is higher 
than the observed value for all r p except at r p = 30/i _1 Mpc. However, as pointed out by J98, the results on scales larger 
than 5/i _1 Mpc are very sensitive to the model for the infall chosen, since the statistical fluctuations become very large. The 
result confirms the well known limit of the CDM model to reproduce the two-point correlation function, to which the PVD is 
connected, since it probes the clustering power on small scales. As stressed several times, the discrepancy of the £(r) on scales 
r p > 5fo _1 Mpc is due to the fact that this model does not have enough power on large scales. Unless the PVD of galaxies is 
biased relative to that of the mass, the CDM model have problems in fitting the observed data. Moreover, using a constant 
bias, it is possible to fit the data only on a limited radius range. Fig. 4 shows that the introduction of a scale-dependent 
bias, due to non-radial motions, produces a notheworthy reduction of the discrepancy of the PVD measured from the Las 
Campanas Redshift Survey with respect to the prediction of the CDM model. The mechanism giving rise to this bias was 
described in the previous section and more widely in DG98, DG99, DG00, D99. Before concluding, it is interesting to note 
that, J98, showed that a scale-dependent bias may explain the discrepancy between the model predictions and observational 
results. In agreement with J98, a velocity bias is not needed to make models compatible with the observed PVD. 
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4 CONCLUSIONS 

In this paper, I have studied the effect of non-radial motions on the PVD in a SCDM model extending the model introduced 
in DG98 and D99. Firstly, I calculated the effect of non-radial motions on the threshold S c and bias. I then compared the large 
scale bias, calculated by means of the model of this paper, with the prediction of Mo & White (1996) model for a spherical 
collapse model and with the result of Sheth et al. (1999) for an ellipsoidal collapse, finding a good agreement with the result 
of this last paper. The model for the bias was used in the SJ97 improved version of CVT to determine the PVD, which was 
compared with the data from DP83, the IRAS value at l/i _1 Mpc of Fisher et al. (1993) and Marzke et al. (1995). As shown 
in Fig. 3, non-radial motions produce a reduction of the values of the PVD, with respect to the prediction of a SCDM model 
that does not take into account this effect. The result is due to the scale-dependent non-local bias induced by non- radial 
motions and to the change of mass distribution in galactic halos. At the light of BB96, the result may be also read as a strong 
dependence of the CVT prediction on the model chosen to describe the mass distribution around galaxies, suggesting that 
the CVT cannot be taken as a direct evidence of a low density universe. 

Finally, I calculated the PVD by means of the redshift distortion in the redshift-space correlation function, £ z (r p ,7r), and 
compared it with the data obtained by J98 from the Las Campanas Redshift Survey. The result confirmed that non-radial 
motions reduce the discrepancy between SCDM model predictions and observations. 
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